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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  can  be  converted  to  SI  (metric)  units  as  follows: 


Multiply 

By 

To  Obtain 

feet 

0.3048 

metres 

feet  per  second  squared 

0.3048 

metres  per  second  squared 

inches 

2.54 

centimetres 

COMBINED  DIFFRACTION  AND  REFLECTION  BY  A  VERTICAL  WEDGE 
PCDFRAC  USER'S  MANUAL 


PART  I:  INTRODUCTTON 


Background 


1 .  In  coastal/ocean  engineering  practice  it  is  often  important  to  be 
able  to  determine  the  wave  heights  near  such  marine  structures  as  jetties, 
breakwaters,  platforms,  and  docks.  Such  information  aids  engineers  in  evaluat¬ 
ing  their  marine  structure  designs,  especially  in  the  areas  of  energy  trans- 
missibility,  sediment  transport,  and  structural  strength.  Wave  heights  in  the 
near  field  have  usua1 ly  been  presented  in  the  form  of  dimensionless  wave 
amplification  factors  which  are  defined  as  the  ratio  of  near  field  wave  height 
to  incident  wave  height. 

2.  Early  studies  (Wiegel  1962)  presented  a  dimensionless  graphical  solu¬ 
tion  for  diffraction  by  a  semi- inf inite  breakwater.  Subsequently,  these 
diagrams  have  been  incorporated  into  every  edition  of  the  Shore  Protection 
Manual  (SPM)  (1984),  and  they  remain  useful  tools  for  preliminary  engineering 
design.  However,  the  diagrams  cannot  account  for  wave  reflection  from  the 
breakwater . 

3.  Chen  (1987)  has  recently  presented  an  analytical  solution  for  the 
total  wave  field  due  to  combined  diffraction  and  reflection  by  a  vertical 
wedge  of  arbitrary  wedge  angle,  as  well  as  the  FORTRAN  program  WEDGE,  which 
calculates  the  amplification  factors  and  the  phases  of  the  waves  in  the  near 
field.  Solutions  for  0-  and  90-deg  wedge  angle,  presented  in  amplification 
factor  diagrams,  are  also  included.  WEDGE  was  designed  to  operate  on  a  con¬ 
ventional  mainframe  computer. 

4  The  FORTRAN  program  PCDFRAC  is  a  version  of  WEDGE  which  is  operable 
on  an  IBM  or  IBM-compatible  personal  computer  (PC).  It  can  he  used  to  calcu¬ 
late  wave  amplification  factors  due  to  combined  diffraction  and  reflection 
near  jettied  harbor  entrances,  quay  walls,  and  other  such  structures.  Jetties 
and  breakwaters  can  be  approximated  by  a  semi-infinite  breakwater  in  which  the 
wedge  angle  is  equal  to  zero.  Corners  of  docks  and  quay  walls  can  be  repre¬ 
sented  by  setting  the  wedge  angle  equal  to  90  deg.  Additionally,  such  natural 
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diffracting  and  reflecting  obstacles  as  rocky  headlands  can  be  approximated  by 
setting  a  particular  value  for  the  wedge  angle. 


5.  The  problem  of  determining  wave  heights  near  marine  structures  has 
been  recognized  for  some  time,  and  it  is  made  more  complex  by  the  combination 
of  water  wave  diffraction  and  reflection  effects  present  near  such  structures. 
Consequently,  a  method  must  be  developed  to  not  only  solve  the  physical  prob¬ 
lem  but  also  accommodate  the  needs  of  various  users  in  a  flexible  and  ef¬ 
fective  manner.  This  report  addresses  this  problem  by  summarizing  an  analyt¬ 
ical  solution  for  the  diffraction  and  reflection  of  monochromatic  incident 
waves  caused  by  a  vertical  wedge.  It  also  presents  PCDFRAC  which  calculates 
the  analytical  solution  to  the  combined  diffraction  and  reflection  problem. 

The  PC  code  makes  the  solution  much  more  available  to  users  than  did  the  main¬ 
frame  version. 

6.  The  objectives  of  this  report  are  to: 

a.  Document  the  program  PCDFRAC. 

b.  Instruct  the  user  (Corps  or  non-Corps)  in  the  proper  applica¬ 
tion  of  the  program  in  field  situations  and  in  the  interpreta¬ 
tion  of  results. 

c.  Instruct  the  user  in  setting  up  PCDFRAC  in  a  personal  computer 
and  to  run  the  code. 

d.  Present  the  user  with  relevant  examples  of  program  execution. 


Report  Organization 


7.  The  report  herein  is  organized  as  follows: 

a.  Part  I  -  introduction  containing  background,  objectives,  and 
program  availability. 

b.  Part  IT  -  brief  theoretical  outline  of  the  basis  of  the  model. 

c.  Part  III  -  details  of  the  model,  including  program  convention 
and  required  equipment  for  proper  program  execution. 

d.  Part  TV  -  procedures  for  installing  the  program  onto  a  personal 
computer . 

e.  Part  V  -  program  execution  and  usage  examples. 
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8.  The  program  PCDFRAC  calculates  wave  heights  due  to  diffraction  and 
reflection  by  a  vertical  wedge  of  arbitrary  wedge  angle  in  constant  water 
depth,  including  the  zero  wedge  angle  (semi-infinite  breakwater)  case.  It 
also  calculates  the  wavelength  from  the  linear  dispersion  relation  by  using  a 
Pade  approximation  (Chen  and  Thompson  1985) .  It  can  be  used  on  an  IBM  or 
IBM-compatible  personal  computer  and  is  available  on  5-1/4-in.*  floppy  disk¬ 
ette  in  Microsoft  FORTRAN  77.  The  diskette  contains  the  source  code 
PCDFRAC. FOR  and  the  executable  file  PCDFRAC.EXE  as  well  as  sample  input  and 
output  files.  It  can  be  obtained  directly  from  the  authors  at  the  following 
address:  Coastal  Oceanography  Branch,  Research  Division,  Coastal  Engineering 

Research  Center,  US  Army  Engineer  Waterways  Experiment  Station,  PO  Box  631, 
Vicksburg,  MS  39180-0631. 


*  A  table  of  factors  for  converting  non-Sl  nntfc  of  measurement  to  ST 
(metric)  units  is  presented  on  page  3. 
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PART  II:  THEORETICAL  BACKGROUND 


9.  This  part  of  the  report  summarizes  the  linear  formulation  and  solu¬ 
tion  of  water  wave  diffraction  and  reflection  by  a  vertical  wedge  of  arbitrary 
wedge  angle,  with  a  semi-infinite  breakwater  considered  a  special  case  in 
which  the  wedge  angle  is  zero.  Water  depth  h*  is  constant,  the  bottom  is 
rigid  and  impermeable,  and  the  monochromatic  incident  waves  of  infinitesimal 
amplitude  come  from  infinity  at  an  angle  a  .  The  cylindrical  coordinate 
system  (r,  6,  z)  is  chosen,  with  z  =  0  representing  the  undisturbed  free 
surface  and  the  positive  z-axis  positioned  vertically  upward.  The  tip  of  the 
wedge  is  chosen  to  be  the  origin  of  the  coordinate  system,  and  the  two  rigid 
walls  of  the  wedge  are  at  6=0  and  6  =  6q  ,  respectively,  as  shown  in 
Figure  1.  Cartesian  coordinates  (x,  y,  z) ,  corresponding  to  the  cylindrical 
coordinates,  are  also  used  on  occasion  and  are  shown  in  the  same  figure.  The 
wedge  angle  is  thus  defined  as  2tt  -  0q  ,  while  the  water  domain  is  defined  as 
0  <  6  i  0  and  0  >  z  >  -  h  . 


*  For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation 
(Appendix  C) . 


10.  For  the  problem  at  hand,  the  velocity  field  for  wave  motion  in  an 
idrji  fluid  may  be  represented  by  the  velocity  potential  <t(r,6,z,t)  and  is 
expressed  as  follows: 


4>(r,0,z,t) 


cosh  k(z  +  h) 
o  cosh  kh 


iu't 

4>  (r , 6 )  e 


( 1 ) 


where 

A  =  -iga  /w 
o  o 

i  =  i-i 

g  =  gravitational  acceleration 
a^  =  incident  wave  amplitude 
w  =  radian  frequency 
k  =  wave  number 

f(r,0)  =  velocity  potential  function  in  the  horizontal  plane 
t  =  time  variable 

The  wave  number  k  must  be  real  and  satisfy  the  following  linear  dispersion 
relation: 


=  gk  tanh  kh  (2) 

11.  Several  properties  of  wave  mechanics  are  dependent  on  the  velocity 
potential  component  ifcfr.O)  .  For  example,  the  free  surface  elevation  n  may 
be  expressed  in  terms  of  f(r,0)  as  follows: 

n  =  a  4(r,P)eiwt  (3) 

o 

The  flow  velocity  u  is  also  dependent  on  4>(r,G)  .  This  quantity  is  impor¬ 

tant  in  the  area  of  sediment  transport  and  may  be  expressed  in  polar  coordi¬ 
nates  as  follows: 


31>  cosh  k(z  +  h)  9(J>  iut 

u  =  —  =  A  - ,  - — n - -  —  e 

r  -r  o  cosh  kii  ar 


(4) 


1  <.*$  cosh  kCz  +  h)  I  3*  lmt 

U ,  =  —  =  - —  =  A  - : — p: -  -  —  e 

0  r  tM  o  cosh  kh  r  3r 


(3) 
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where  the  subscripts  r  and  6  refer  to  the  r-  and  6-directions .  The  abso¬ 
lute  values  of  u  and  un  are  the  maximum  flow  velocities  for  each  com- 
r  0 

ponent  direction,  while  the  phases  of  u  and  u^  contain  the  wave  phase 
information.  It  is  clear,  then,  that  the  solution  to  Equation  1,  as  well  as 
to  Equations  3  through  5,  lies  in  finding  ij>(r,0)  • 

12.  For  an  incident  plane  wave  train  coming  from  the  a  direction  (as 
shown  in  Figure  1)  where  the  free  surface  elevation  of  the  incident  wave  may 
be  described  by: 


i(kr  cos  a  +  iot) 

n .  =  a  e 
l  o 

the  analytical  solution  for  a  wave  field  by  a  vercical  wedge  of  arbitrary 
wedge  angle  (based  on  the  linearized  wave  theory)  may  be  written  as  follows 
(Chen  1987): 


<p(r,9) 


J  (kr) 
o 


+  2 


2 

n=l 


inTr/2v  T 
e  J 


n/ v 


na  n0 

cos  —  cos  — 
V  V 


(7) 


where 

V  =  0  /  7T  (8) 

o 

J  =  zeroth  order  Bessel  function  of  the  first  kind 
o 

J  ,  =  n/v  order  Bessel  function  of  the  first  kind 

n/v 

The  semi-infinite  breakwater  is  a  special  case  of  the  diffraction  and  reflec¬ 
tion  problem  where  the  wedge  angle  is  equal  to  zero  and  v  =  2  .  The  solution 
of  Equation  1  for  this  case  is 


<Kr,0) 


T  n  -I  ,  O  V  in'M  Jn  na  n. 

J  (kr)  +  2  >  e  tt—  cos  —  cos  — 

o  2  2 

n=  I 


f  93 


The  velocity  potential  function  <j>(r,0)  in  Equations  7  and  9  is  a  complex 
function  and  may  be  expressed  as 


<j>  =  I  <j>  |  e 


where 


£l;n(Oj“  +  [Re(<£)J“  =  amplitude  of 


6  =  tari_1  [iffry] =  phase  of  * 


and  where 


Tm  <(>  =  imaginary  part  of  $ 

Re  if  =  real  part  of  <f> 

Substituting  Equation  10  into  Equation  3,  the  following  is  obtained: 


1(6+ cot) 


n  =  a^  |  <p  [  e 


This  expression  represents  the  actual  water  surface  elevation  at  a  point  in 
the  water  domain  bounded  by  a  vertical  wedge  of  arbitrary  wedge  angle.  Since 
tne  incident  wave  train  is  expressed  in  Equation  6,  the  normalized  water  sur¬ 
face  elevation  in  the  near  field  may  be  expressed  as 

i(6-kr  cos  a) 


where  is  the  incident  free  surface  elevation.  It  is  clear  that  the  term 


i(B-kr  cos  a) 


is  a  factor  which  modifies  the  incident  wave  elevation  to 


account  for  reflection  and  diffraction  effects;  thus,  the  amplitude  of  the 
normalized  surface  elevations  may  be  expressed  as  the  following  wave  amplitude 
factor: 


The  phase  of  Equation  14  is  the  following  phase  difference  between  incident 
and  amplified  waves: 
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(16) 


Phase  of  —  =  6  -  kr  cos  a 
ni 

13.  Output  of  the  program  PCDFRAC  is  comprised  of  |  <J>  |  and  8  .  The 
amplification  factor  |  |  would  then  be  multiplied  by  the  incident  wave 
height  to  obtain  the  actual  wave  height.  The  phase  of  the  amplified  wave  8 
is  a  quantity  not  usually  required  in  engineering  practice;  however,  it  may  be 
useful  on  some  occasions. 
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PART  III:  PCD FRA C  OPERATIONAL  DETAILS 


Program  Characteristics 


14.  The  program  PCDFRAC  has  been  developed  specifically  for  the  problem 
of  water  wave  diffraction  and  reflection  from  a  vertical  wedge  of  arbitrary 
wedge  angle.  It  has  been  designed  for  operation  on  an  IBM  or  IBM-compatible 
PC.  A  thin,  rigid,  impermeable  semi-infinite  breakwater  is  a  special  case  in 
which  the  wedge  angle  is  equal  to  zero.  The  program  language  is  FORTRAN  77. 
The  source  code,  PCDFRAC. FOR,  requires  114  Kb  of  space,  while  the  executable 
file,  PCDFRAC.EXE,  uses  74  Kb. 


Input  and  Output 

15.  The  program  is  divided  into  several  parts.  The  first  part  is  the 
main  program  PCDFRAC,  with  which  the  user  works  directly.  It  calls  the  solv¬ 
ing  subroutines  and  prints  the  answers  in  usable  form.  It  also  calculates 
wavelength  from  the  input  wave  period  by  the  Pade  approximation  to  the  linear 
dispersion  relation  (Chen  and  Thompson  1985).  The  subroutine  PCWEDGES  solves 
the  actual  problem  and  uses  the  mathematical  subroutines  BESJ,  JAIRY,  and 
GAMLN  which  are  borrowed  from  the  Naval  Surface  Weapons  Center  (NSCW)  Library 
of  Mathematics  Subroutines  (Morris  1984).  Onlv  the  main  program  PCDFRAC  and 
the  subroutine  PCWEDGES  are  listed  in  Appendix  B.  The  main  code  of  PCDFRAC  is 
similar  to  that  of  WEDGE  (See  Chen  1987). 

16.  Inputs  into  the  program  are  as  follows: 

a.  Constant  water  depth  D  ,  in  feet. 

b.  Period  T  of  the  incident  wave,  in  seconds. 

£.  Approach  angle  WAVEA  of  the  incident  wave,  in  degrees 
counter-clockwise  from  the  positive  X-axis. 

d.  Wedge  angle  WEDGEA  of  the  obstacle,  in  degrees.  A  default  of 
the  program  allows  the  semi-infinite  breakwater  case  (WEDGEA 
=  0.0  deg)  to  be  run. 

e.  Location  X,Y  of  the  desired  calculation,  in  feet  from  the  tip 
of  the  wedge. 

Figure  2  shows  the  preceding  input  conventions. 


I? 


17.  The  program  outputs  are: 

a.  Wavelength  L  ,  In  feet. 

b.  Wave  amplification  factor  ]<J>|  as  given  in  Equation  11,  which 
is  the  ratio  of  wave  height  to  incident  wave  height. 

£.  Wave  phase  6  as  given  in  Equation  12,  in  radians. 

18.  The  program  may  be  run  either  interactively  (i.e.,  working  directly 

with  the  user),  or  as  a  batch  job  (where  the  data  is  read  from  a  file).  These 
features  will  be  discussed  in  further  detail  in  Part  V. 

19.  Equations  7  and  9  showed  a  summation  of  an  infinite  number  of 
terms.  This  summation  has  been  accommodated  in  the  program  by  carrying  the 

summation  out  to  a  term  followed  by  eight  successive  terms  in  which  the  abso- 

—8 

lute  value  of  the  Bessel  function  is  10  or  less.  If  the  value  of  the  solu- 

—8 

tion  is  of  the  order  of  one,  this  corresponds  to  a  truncation  error  of  10  or 
less . 

20.  For  values  of  R/L  (radial  distance  to  point  of  calculation 
divided  by  wavelength)  greater  than  16,  the  statement  "NO.  OF  TERMS  FOR  SUM¬ 
MATION  IS  INSUFFICIENT"  may  appear.  This  statement  indicates  that  the  number 

_8 

of  terms  in  the  summation  is  insufficient  to  ensure  the  10  truncation  error 
and  that  the  NN  value  in  the  PARAMETER  statement  in  the  subroutine  PCWEDGES 
needs  to  be  greater  than  the  present  value  of  200  if  accuracy  is  to  be 
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maintained.  In  this  case,  the  user  should  change  NN  from  200  to  a  larger 
value . 

21.  A  program  of  this  length  would  require  a  certain  amount  of  run  time 
on  a  PC.  Tables  A1  and  A2  (located  in  Appendix  A)  give  run  times  (in  seconds) 
for  various  values  of  the  truncation  error  and  R/L  ratio.  The  times  in¬ 
dicated  are  the  lengths  of  time  the  computer  would  require  to  process  one  set 
of  (x,y)  coordinates  for  one  incident  wave  angle  and  one  incident  wave  period 
(on  the  order  of  seconds) .  It  follows  that  the  greater  the  truncation  error 
allowed,  the  less  time  the  computer  would  need  to  process  one  set  of  coordi¬ 
nates.  It  also  follows  that  the  greater  the  R/L  ratio,  the  longer  the  com¬ 
puter  would  take  to  finish  processing  one  set  of  coordinates,  as  more  terms  in 
the  summation  are  required  to  ensure  proper  accuracy.  Table  AI  is  a  collec¬ 
tion  of  run  times  using  an  IBM  AT-class  computer,  while  Table  A2  is  comprised 
of  run  times  on  an  IBM  XT  computer. 


Required  Equipment 


22.  The  source  code  PCDFRAC.FOR  and  the  executable  file  1DFRAC.EXE  are 
contained  on  a  floppy  diskette  which  may  be  obtained  upon  user  request.  Addi¬ 
tionally,  the  following  equipment  is  either  required  or  recommended: 

a.  An  IBM  or  IBM-compatible  personal  computer  with  a  hard  disk  and 
disk  operating  system  (DOS) . 

b.  A  FORTRAN  compiler  (installed  in  the  computer) .  Major 
differences  between  various  brands  of  compilers  are  not  ex¬ 
pected;  however,  the  compiler  must  accept  the  Fortran  77 
language . 

c.  One  5-1/4  in.  disk  drive. 

d.  Printer  (optional) . 

e.  Text  editor. 
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PART  IV:  START-UP  PROCEDURES 


Copying  Program  onto  Hard  Disk 

23.  As  mentioned  earlier  in  this  report,  an  IBM  or  IBM-compatible  per¬ 
sonal  computer  with  DOS  is  required  to  use  this  program.  Although  Microsoft 
DOS  Version  2.10  was  used  for  this  copying  procedure,  major  deviations  among 
various  brands  of  DOS  are  not  expected.  To  perform  this  procedure  the  user 
must : 

a.  Insert  the  floppy  disk  into  a  particular  disk  drive,  and  switch 
the  operating  drive  of  the  computer  to  the  drive  and  directory 
containing  the  Fortran  compiler.  For  Microsoft  DOS  this  is 
done  by  typing 

C:  CD /FORTRAN  <CR> 


if  the  Fortran  compiler  were  located  in  drive  C  and  in 
directory  FORTRAN.  (The  symbol  <CR>  denotes  carriage  return.) 

b.  Copy  the  program  PCBFRAC.FOR  onto  the  hard  disk.  For 
Microsoft  DOS  Version  2.10  this  is  done  by  typing 

COPY  A:  PCDFRAC .FOR  <CR> 


if  the  disk  were  in  drive  A.  Now  the  program  PCDFRAC. FOR  is 
copied  onto  the  same  disk  and  into  the  same  directory  as  the 
FORTRAN  compiler. 


Compilation  of  Program 

24.  Compilation  of  a  program  translates  the  language  in  which  the  pro¬ 
gram  is  written  (in  this  case,  FORTRAN)  into  the  internal  machine  language  of 
the  computer.  The  procedure  here  is  that  of  the  Ryan-McFarland  Version  1.00 
Fortran  compiler;  however,  major  discrepancies  among  various  compilers  are  not 
expected.  To  compile  the  program,  the  user  should  type 

PROFORT  PCDFRAC/E  <CR> 


and  compilation  will  begin.  With  the  Ryan-McFarland  Version  1 .00  compiler. 
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the  /E  option  added  to  the  end  of  the  file  name  will  list  any  syntactical 
errors  as  they  occur.  There  should  be  no  errors. 

25.  Compiling  the  program  on  an  IBM  XT  or  XT-compatible  computer  takes 
about  12  min,  primarily  because  of  the  length  of  the  program.  On  an  IBM  AT  or 
AT-compatible  personal  computer,  compilation  requires  just  over  4  min.  When 
the  compilation  process  is  complete,  the  following  message: 

Compilation  Complete:  0  Errors,  0  Warnings 

should  appear  on  the  screen.  An  object  file,  PCDFRAC.OBJ,  has  been  created  as 
a  result. 

Linkage  of  Program 

26.  Linkage  of  the  program  with  any  external  subroutines  is  done  next. 
This  is  a  necessary  step  with  all  programs,  regardless  whether  any  external 
subroutines  are  specified  within  the  program.  Linkage  is  done  by  typing 

LINK  PCDFRAC  <CR> 

The  following  will  then  print  to  the  screen,  one  at  a  time  (each  to  be 
answered  with  a  carriage  return): 


RUN  FILE 

[PCDFRAC.EXE] 

<CR> 

LIST  FILE 

[NUL.MAP] 

<CR> 

LIBRARIES 

[.LIB] 

<CR> 

The  above  options  actually  refer  to  specific  options  in  the  linking  process. 
These  options  are  not  directly  involved  with  PCDFRAC  and  are  not  required  for 
program  execution.  Linkage  takes  about  30-40  sec,  regardless  of  the  type  of 
computer.  After  linkage,  the  usual  computer  prompt  will  appear,  and  an  execu¬ 
table  file,  PCDFRAC.EXE,  will  have  been  created. 


•: 


i 
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Execution  of  Program 

27.  The  executable  file,  PCDFRAC.EXE,  actually  runs  the  program.  If 
the  user  wished  to  use  PCDFRAC  without  any  Intent  to  modify  It,  this  file 
would  be  used  exclusively.  Additionally,  no  FORTRAN  compiler  is  required  to 
use  PCDFRAC.EXE.  To  perform  any  modifications,  however,  the  file  PCDFRAC. FOR 
would  be  needed.  If  the  program  were  modified  in  any  way,  it  would  need  to  be 
recompiled  and  relinked  after  the  changes  were  completed.  To  execute  the  pro¬ 
gram  the  user  should  type 

PCDFRAC  <CR> 


and  the  program  will  begin  running. 


PART  V:  PROGRAM  EXECUTION 

28.  The  program  will  offer  the  option  of  either  operating  directly  with 
the  user  (interactive  operation)  or  of  reading  and  processing  the  data 
directly  from  an  existing  file  (batch  operation) . 

Interactive  Option 

29.  The  interactive  option  is  recommended  for  users  with  a  small  amount 
of  data.  It  is  sequenced  as  follows: 

a.  The  program  will  first  ask  whether  a  written  output  is  desired. 

If  so,  the  user  will  input  the  file  name  to  be  used.  The 
program  contains  a  feature  which  will  alert  the  user  if  the 
file  name  specified  matches  that  of  an  existing  file.  If  this 
occurs,  the  program  will  offer  the  user  the  option  of  either 
entering  a  new  file  name  or  of  overwriting  the  existing  file. 
(Note:  Limited  to  a  maximum  of  an  eight-character  string,  the 

file  name  can  have  fewer  characters.  Some  examples  of 
eight-character  strings  include:  FAST.DAT,  DIFRACT1 ,  and 
WAVES- 10.) 

b.  The  program  will  then  ask  if  the  input  data  are  stored  in  a 
file.  Because  the  interactive  option  is  to  be  used  here,  the 
answer  is  "no". 

c.  The  program  will  then  ask  if  the  user  would  require  a 
particular  value  for  the  wedge  angle  or  the  semi- inf inite 
breakwater  case.  If  a  value  for  the  wedge  angle  is  needed,  it 
would  be  input  (in  degrees)  during  this  step. 

d.  The  program  will  then  request  the  following: 

Water  depth,  in  feet. 

Incident  wave  period,  in  seconds. 

Incident  wave  angle,  in  degrees  counter-clockwise  from  the 
positive  x-axis. 

e.  Next  the  program  will  ask  for  the  number  of  points  to  be  input 
on  the  present  run. 

f_.  It  will  then  ask  for  the  (x,y)  coordinates  (in  feet)  of  the 
desired  points  of  calculation,  one  at  a  time. 

g.  After  all  the  points  have  been  entered,  it  will  then  give  the 
information  in  tabular  form,  including  the  input  information. 

If  the  output  were  also  written  to  a  file,  the  user  would  be 
reminded  of  the  file  name. 


18 


30.  In  the  example  problem  shown  in  Figure  3,  it  is  desired  to  cal¬ 
culate  the  wave  amplification  factors  at  the  locations  shown,  with  the 
obstacle  being  a  semi-infinite  breakwater  (wedge  angle  =  0.00).  The  inter¬ 
active  run  for  the  example  in  Figure  3  is  shown  in  Figure  4  with  user  input 
underlined.  The  resulting  output  is  then  sent  to  both  the  screen  and  any 
output  file  specified.  It  appears  in  Figure  5. 


v 


BREAKWATER 

(WEDGE  ANGLE  =0.0  DEG.) 


•  •  * 

(-50,-25)  (0,-25)  (50,-25) 

Figure  3.  Example  problem  for  semi-infinite  breakwater 


Explanation  of  Output 


31.  In  the  interactive  output  shown  in  Figure  5,  the  first  three 
columns,  as  well  as  the  fifth,  sixth,  and  seventh  columns,  consist  of  user  in¬ 
put  which  require  no  further  explanation.  The  fourth  column  is  the  calculated 
wavelength.  The  last  two  columns  comprise  the  primary  solution,  as  given  in 
Equations  11  and  12.  AMP  is  the  amplification  factor  for  the  incident  wave. 
This  number  would  be  multiplied  by  the  incident  wave  height  to  obtain  the 
actual  wave  height  at  that  particular  (X,Y)  location.  PHA  is  the  wave  phase 
(in  radians).  This  information  is  not  usually  required  in  most  engineering 
applications  but  may  be  of  use  at  times. 
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*************** ********************************************* ********** 


*  * 

*  "PCDFRAC"  -  A  PROGRAM  WHICH  CALCULATES  WAVE  AMPLIFICATION  FACTORS  * 

*  BASED  ON  THE  COMBINED  DIFFRACTION  AND  REFLECTION  OF  WAVES  BY  A  * 

*  VERTICAL  WEDGE  OF  ARBITRARY  WEDGE  ANGLE.  THE  PROGRAM  ALSO  CALCU-  * 

*  LATES  THE  PHASE  OF  THE  AMPLIFIED  WAVE.  THE  PROGRAM  WILL  DEFAULT  TO  * 

*  ZERO  WEDGE  ANGLE  ( SEMI - INFINITE  BREAKWATER  CASE )  IF  NO  WEDGE  ANGLE  * 

*  IS  SPECIFIED.  * 

*  * 


********************************************************************** 


DO  YOU  WANT  OUTPUT  IN  A  FILE  ?  (0="NO",  1="YESM) 


1 


ENTER  NAME  OF  OUTPUT  FILE  =  TEST  I . OP 

IS  INPUT  DATA  LOCATED  IN  A  FILE  ?  ( 0= 'NO'  , 1  =  " YES" )  =  0_  • 

OPTION  TO  SET  WEDGE  ANGLE  (DEFAULT  VALUE  -  0.00  DEG.) 

SET  WEDGE  ANGLE  ?  (  0  -  "NO "  ,  1  = "  YES"  )  =  _0_ 

SET  INITIAL  VALUES 

INPUT  WATER  DEPTH  (FT.  )  =  _20 

INPUT  INCIDENT  WAVE  PERIOD  (SEC.)  -  8_ 

INPUT  INCIDENT  WAVE  ANGLE 
i DEG .  CCW  FROM  POSITIVE  X  -  AXIS)  -  135 


HOW  MANY  POINTS  DO  YOU  WISH  TO  ENTER  (UP  TO  99)?  -  _9 

AT  EACH  PROMPT,  DO  THE  FOLLOWING  : 

1)  TYPE  IN  X  -  COORDINATE 
~ )  TYPE  A  COMMA 

3)  TYPE  IN  Y  -  COORDINATE 

4)  IRES 3  ENTER" 


X.  Y) 

COORDINATES 

F  OH 

POINT 

NO. 

1 

(FT.  ■ 

-  -50.25 

V  V 

COORDINATES 

F:  >R 

I  "'I NT 

NO. 

- 

( FT .  ' 

—  o  _  3  k 

X.  Y) 

COORDINATES 

For 

POINT 

NO. 

3 

(  FT.  ) 

-  50,25 

X ,  Y 

COORDINATES 

FOR 

PJINT 

NO. 

4 

(  FT.  ) 

<LT- 

O 

o 

X,  Y) 

COORD  I  NATES 

FOR 

I'OINT 

NO 

!;> 

(FT.  ) 

-  0 , 0 

X .  Y ) 

COORDINATES 

F  OR 

POINT 

NO. 

6 

(  FT  .  j 

-  50,0 

X .  i  ) 

COORDINATES 

Fok 

POINT 

NO . 

i 

(  FT .  i 

-  -50,-2! 

X ,  Y  j 

COORDINATES. 

FOR 

POINT 

NO. 

8 

(  FT.  ) 

-  0.-25 

X  .  Y  i 

COORDINATES 

FOR 

POINT 

NO . 

9 

■  FT .  < 

-  50 , -25 

Figure  4.  Interactive  execution  for  problem  in  Figure  3 
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SEMI- INFINITE  BREAKWATER  (WEDGE  ANGLE  =0.00  DEGREES) 


DEPTH 

PERIOD 

WAV . ANGLE 

WAVLNG . WDG . ANG . 

X 

Y 

AMP 

PHA 

(FT) 

(SEC) 

(DEG) 

(  FT )  ( DEG  ) 

(FT) 

(FT) 

-  -  - 

(RAD) 

20.00 

8.00 

135.00 

190.00  0.00 

-50 . 00 

25.00 

0.98 

1 . 77 

0.00 

25.00 

1 . 04 

0.45 

50 . 00 

25.00 

1 . 44 

-0 . 81 

-50 . 00 

0.00 

1 . 00 

1  .  17 

0 . 00 

0 . 00 

1 . 00 

0.00 

50.00 

0 . 00 

1  .  69 

-0.90 

-50 . 00 

-25 . 00 

1.03 

0 . 59 

0 . 00 

-25.00 

0.86 

-0 . 55 

50.00 

-25.00 

0 . 60 

-2.18 

YOUR  OUTPUT  IS  IN  TEST  I. OP 
Execution  terminated  :  0 


C> 

Figure  5.  Solution  to  problem  in  Figure  3 

32.  Figure  6  is  an  example  problem  in  which  the  diffracting  obstacle  is 
a  dock  with  a  90-deg  wedge  angle.  It  is  desired  to  calculate  the  wave  ampli¬ 
fication  factors  at  the  locations  shown.  The  interactive  run  and  the  output 
are  both  similar  to  those  of  the  example  shown  in  Figure  3  and  require  no 


y 
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further  explanation.  The  interactive  run  is  shown  in  Figure  7,  while  the  cor¬ 
responding  output  is  shown  in  Figure  8. 

Batch  Option 

33.  The  batch  operation  is  recommended  for  calculation  of  a  large 
number  of  points  or  if  the  user  requires  calculations  with  several  values  of 
wave  period,  incident  wave  angle,  and/or  wedge  angle.  The  batch  operation  is 
sequenced  as  follows: 

a.  The  user  will  first  create  the  file.  An  example  appears  in  a 
following  subsection.  (Note:  The  file  name  must  be  kept  to  a 
maximum  of  eight  characters,  as  before.) 

b.  The  program  will  ask  if  a  written  output  file  is  desired.  If 
it  is,  the  user  will  enter  the  name  of  the  output  file.  As 
before,  the  program  will  alert  the  user  if  the  chosen  output 
file  name  matches  that  of  an  existing  file.  If  so,  the  user 
will  have  tne  option  of  either  overwriting  the  file,  or  er.ter- 
a  new  file  name. 

c.  it  will  then  ask  if  the  data  are  located  in  a  file.  Because 
batch  operation  is  desired,  the  answer  input  bv  the  user  would 

i  I!  M 

be  yes  . 

d.  It  will  then  request  the  input  file  name. 

e.  As  it  calculates,  it  will  print  the  results  to  the  screen  in 
tabular  form.  if  a  written  file  was  requested,  it  will  write 
to  the  file  in  the  same  format. 

_f .  When  completed,  it  will  ask  the  user  it  nnv  additional  files 
need  to  be  processed.  If  so,  the  user  will  input  the  name  of 
t  he  f  i  !  e  . 

£.  The  entire  sequence  is  then  repeated  until  there  ate  no  more 
input  files  left  to  process. 

34.  The  input  f i 1 e  for  the  batch  job  requires  a  specific  formal  and  can 
be  created  with  any  compatible  text  editor.  The  input  file  in  Figure  4 
contains  data  lor  the  example  problems  in  Figures  3  and  6.  There  is  a  limit 
of  3<10  values  of  X  and  Y  which  can  be  processed  in  one  batch  job.  However, 

a  greater  number  of  values  can  be  processed  bv  changing  the  value  of  300  in 
the  DIMENSION  statement  in  program  PCPFRAC . FOR  to  a  higher  value.  To  create 
an  input  file,  the  following  info) mat  ion  is  used: 

a .  Line  1  -  P,T,WAV FA.WF.DO -FA 

The  water  depth,  incident  wave  period,  incident  wave  angle  and 
wedge  angle  are  input  here.  D  appears  in  columns  l-’,0,  \  in 
columns  Ji-20,  WAVF.A  in  columns  21-30.  and  WF.DGFA  in 


********************************************************************** 
*  * 

*  "PCDFRAC"  -  A  PROGRAM  WHICH  CALCULATES  WAVE  AMPLIFICATION  FACTORS  *■ 

*  BASED  ON  THE  COMBINED  DIFFRACTION  AND  REFLECTION  OF  WAVES  BY  A  * 

*  VERTICAL  WEDGE  OF  ARBITRARY  WEDGE  ANGLE.  THE  PROGRAM  ALSO  CALCU-  * 

*  LATES  THE  PHASE  OF  THE  AMPLIFIED  WAVE.  THE  PROGRAM  WILL  DEFAULT  TO  * 

*  ZERO  WEDGE  ANGLE  . SEMI  -  INFINITE  BREAKWATER  CASE)  IF  NO  WEDGE  ANGLE  * 

*  IS  SPECIFIED.  * 

*  * 
********************************************************************** 


DO  YOU  WANT  OUTPUT  IN  A  FILE  ?  (0="NO",  1=  'YES" )  -  1 


ENTER  NAME  OF  OUTPUT  FILE  -  TESTI.OP 


IS  INPUT  DATA  LOCATED  IN  A  FILE  ?  (0="NO“ , 1="YES" )  -  0. 


OPTION  TO  SET  WEDGE  ANGLE  ( DEFAULT  VALUE  =  0.00  DEG.  i 
SET  WEDGE  ANGLE  7  i 0=  NO" . 1= "YES" )  =  1 
ENTER  WEDGE  ANGLE  (DEG.)  =  90 


SET  INITIAL  VALUES 


INPUT  WATER  DEPTH  (FT.)  =  20 
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k  SEC .  ) 
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TO  99  i?  =  _3_ 

AT  EAC 

H  PROMPT 

00 
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0 PLOWING 
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x  -  • 
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NATE 
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TYPE  A  0 
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4  ) 
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COORD: NA 
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INT 
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(FT. 

=  -50.25 

X .  Y 

bbUkLi  x  NA 

TEb  X 

OR  l  0 
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NO 
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k  FT . 

=  0.25 

X ,  V 

’OORT INA 

TES  ? 
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INT 
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;  pr 

-  f,  0  t  2  5 

X .  V 

COORDI NA 

TES  F 
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INT 
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(  FT 

=  -50,0 

X ,  Y 

'uORI-'  1 1 J  A 

TEb  r 
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INI 
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..  FT. 

=  0.0 
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Figure  7.  Interactive  execution  for  problem  in  Figure  (3 
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OBSTACLE  IS  A  VERTICAL  WEDGE  OF  ANGLE  90.00  DEGREES 


DEPTH 

PERIOD 

WAV. ANGLE 

WAVLNG . 

WDG . ANG . 

X 

Y 

AMP 

PHA 

(FT) 

( SEC ) 

i  DEG  ) 

(FT) 

(DEG) 

1  FT  ) 

..FT) 

_ 

(RAD) 

20.00 

8 . 00 

135 . 00 

190.00 

90 . 00 

-50 . 00 

25.00 

0.30 

1 . 94 

0.00 

25 . 00 

1.01 

0 . 30 

50.00 

2  5.00 

1.43 

-0.87 

-50.00 

0 . 00 

0 . 75 

: .  18 

0.00 

0.00 

1  .  33 

0 . 00 

50.00 

0.00 

1  .  71 

-0 . 95 

-50.00 

-  25 . U0 

0.75 

0  .  35 

0 . 00 

-25 . 00 

1  .  53 

-0.41 

YOUR  OUTPUT  IS  IN  TEST I. OP 
Execution  terminated  :  0 


C> 


Figure  8.  Solution  to  problem  in  Figure  6 


20.00 

3.00 

135.00 

90 . 00 

0.00 

-25.00 

-50.00 

-25.00 

50.00 

0 . 00 

0 . 00 

0.00 

-50 . 00 

0.00 

50.00 

25.00 

0.00 

25.00 

50 . 00 

25 . 00 

1 

20 . 00 

8.00 

135.00 

50.00 

-25.00 

0.00 

-25 . 00 

-50.00 

-25.00 

50 . 00 

0 . 00 

0  .  U0 

0 . 00 

-50.00 

0.00 

50.  U0 

25.00 

0.00 

25.00 

-50.00 

25.00 

1 

Figure  9.  Sample  input  data  file  for  Figures  3  and  6 

columns  31-40.  If  the  semi-infinite  breakwater  rase  is  to  he 
run,  the  input  field  for  WEDGEA  can  be  left  blank. 

b.  Line  2  -  X , Y , END 

The  location  of  the  calculation  is  input  on  each  line.  X  is 
input  in  columns  1-10,  Y  in  columns  11-20,  and  END  in 
column  50.  END  is  an  end-of-file  marker  alerting  the  computer 
that  that  particular  line  is  the  last  line  of  X,Y  data  for  the 
particular  values  of  D,  T,  WAVEA,  and  WEDGF.A.  To  denote  the 
end  of  X,Y  data,  the  user  should  enter  1  for  END.  Otherwise, 
the  input  field  may  be  left  blank.  Successive  sets  of  D,  T, 
WAVEA,  and  WEDGF.A  data  can  then  be  entered,  along  with 
their  correspond ing  sets  of  X,Y  data. 


35.  The  execution  of  the  data  file  in  Figure  9  is  shown  in  Figure  10 
with  the  user  input  underlined.  Figure  11  shows  the  resultant  output. 


********************************************************************** 


*  * 

*  PCDFRAC  -  A  PROGRAM  WHICH  CALCULATES  WAVE  AMPLIFICATION  FACTORS  * 

*  BASED  ON  THE  COMBINED  DIFFRACTION  AND  REFLECTION  OF  WAVES  BY  A  * 

*  VERTICAL  WEDGE  OF  ARBITRARY  WEDGE  ANGLE.  THE  PROGRAM  ALSO  CALC V-  * 

*  LATES  THE  PHASE  OF  THE  AMPLIFIED  WAVE.  THE  PROGRAM  WILL  DEFAULT  TO  * 

*  ZERO  WEDGE  ANGLE  ( SEMI - INFINITE  BREAKWATER  CASE)  IF  NO  WEDGE  ANGLE  * 

*  IS  SPECIFIED.  * 

*  * 


************************  ****************** **************  ***********  *** 


DO  YOU  WANT  OUTPUT  IN  A  FILE 

?  (0=  "NO" , 

1 = "YES" ) 

-  J_ 

ENTER  NAME  OF  OUTPUT  FILE  = 

TESTB.OP 

IS  INPUT  DATA  LOCATED  IN  A  FILE  ?  (0-"NO" , 1="YES" )  -  _L 
ENTER  NAME  OF  INPUT  FILE  =  TESTB.IP 


Figure  10.  User  input  for  execution  of  file  in  Figure  9 
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INPUT  FILE  -  TESTB. IP 


DEPTH 

PERIOD 

WAV. 

ANGLE 

WAVLNG . 

WDG 

.  ANG 

X 

Y 

AMP 

PHA 

(FT) 

(SEC) 

(DEG) 

(FT) 

(DEG) 

l 

FT) 

(. 

FT) 

_ 

_  __ 

1 

RAD) 

20. 

.00 

3.00 

135 

.00 

190 

.  00 

90 

.00 

0 

.  00 

-25 

.  00 

1 

.  53 

-0 

.41 

-50 

.  00 

-25 

.00 

0 

.  75 

0 

.  35 

50 

.00 

0 

.00 

1 

.71 

-0 

.  95 

0 

.00 

0 

.00 

1 

.  33 

0 

.00 

-50 

.00 

0 

.  00 

0 

.  75 

i 

.  18 

50 

.  00 

25 

.00 

1 

.43 

-0 

.  87 

0 

.00 

25 

.00 

1 

.01 

0 

.  30 

-50 

.00 

25 

.00 

0 

.  90 

1 

.  94 

20. 

.00 

8.00 

135 

.00 

190. 

00 

0. 

00 

50 

.00 

-25 

.00 

0 

.  60 

-2 

.  18 

0 

.00 

-25 

.00 

0 

.  86 

-0 

.  55 

-50 

.00 

-25 

.00 

1 

.  03 

0 

.  59 

50 

.00 

0 

.00 

1 

.  69 

-0 

.  90 

0 

.00 

0 

.00 

1 

.  00 

0 

.00 

-50 

.00 

0 

.00 

1 

.  00 

1 

.  17 

50 

.  00 

25 

.00 

1 

.  44 

-0 

.  81 

0 

.00 

25 

.00 

1 

.0-, 

0 

.  45 

-50 

.00 

25 

.00 

0 

.  98 

1 

.  77 

YOUR  OUTPUT  IS  IN  TESTB. OP 

DO  YOU  HAVE  ANOTHER  INPUT  FILE  ?  (0="NO  " , 1-  "YES" )  = _0 
Execution  terminated  '•  0 


C> 


Figure  11.  Output  for  execution  of  sample  input  file  in  Figure  9 
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PART  VI :  SUMMARY 


36.  A  brief  theoretical  outline  of  wave  diffraction  and  reflection  by  a 
vertical  wedge  of  arbitrary  wedge  angle  has  been  presented.  A  personal 
computer  code,  PCDFRAC,  which  is  a  PC  version  of  WEDGE  (Chen  1987),  has  been 
implemented  for  calculating  the  wave  amplification  factors  necessary  to  obtain 
wave  heights  in  the  near  field  of  the  vertical  wedge  due  to  diffraction  and 
reflection . 

37.  The  necessary  procedures  for  setting  up  PCDFRAC  on  a  personal 
computer  have  been  documented  in  detail.  Input  and  output  files  have  also 
been  described  in  detail.  Possible  applications  of  the  program  in  field 
situations  have  been  discussed  and  the  execution  of  the  program  outlined  in 
sequential  order.  Instructive  examples  of  the  execution  of  PCDFRAC  have  also 
been  presented. 
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APPENDIX  Ar  RUN  TIMETABLES 


Table  A1 

Run  Time*  Versus  Summation  Truncation  Error 
IBM  AT-Class 


Truncation  Error 

R/L  Ratio 

E-08 

E-06 

E-04 

1.0 

0.82 

0.70 

0.59 

2.0 

2.42 

2.26 

1.98 

3.0 

2.56 

2.28 

2.06 

4.0 

2.72 

2.31 

2.10 

5.0 

3.29 

2.85 

2.46 

6.0 

3.41 

3.06 

2.71 

7.0 

3.91 

3.46 

3.08 

8.0 

4.29 

3.88 

3.37 

9.0 

4.84 

4.28 

3.83 

10.0 

5.20 

4.70 

4.10 

*  In  seconds. 

Run  Time* 

Table  A2 

Versus  Summation 

Truncation  Error 

IBM 

XT 

Truncation  Error 

R/L  Ratio 

E-08 

E-06 

E-04 

1.0 

2.15 

1.65 

1 .4Q 

2.0 

4.20 

4.11 

3.03 

3.0 

4.85 

4.76 

4.29 

4.0 

5.65 

5.08 

4.32 

5.0 

5.85 

5.10 

4.38 

6.0 

6.04 

5.65 

4.50 

7.0 

7.04 

6.04 

5.50 

8.0 

7.25 

6.60 

5.97 

9.0 

8.14 

6.87 

6.10 

10.0 

8.95 

7.66 

7.17 

*  In  seconds. 


« 
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APPENDIX  B:  PARTIAL  LISTING  OF  PROGRAM  PCD FRA C 


PROGRAM  PCDFRAC 

C  **************************************************** 

C  *This  program  calculates  wave  amplification  factors* 


C  *for  the  combined  diffraction  and  reflection  of  * 
C  ^monochromatic  incident  waves  of  infinitesimal  * 
C  ^amplitude  coming  from  infinity  by  a  vertical  wedge* 
C  *of  arbitrary  wedge  angle.  It  also  calculates  the  * 
C  *phase  of  the  amplified  wave.  The  program  will  * 
C  *default  to  a  zero  wedge  angle  i semi-infinite  * 
C  *breakwater)  if  no  wedge  angle  is  specified.  * 


C  **************************************************** 

C 

C  Initialize  and  dimension  variables 

C 

COMPLEX  F 
REAL  L.PI.PI2 
INTEGER  END 
CHARACTER* 8  OUTFILE 
CHARACTER*8  INFILE 
LOGICAL  EXST 

DIMENSION  X( 300 ) , Y( 300 ) , FABS( 300 ) , FPHA( 300 ) 

PI=3. 14159 
P12-2 . *PI 
WEDGEA=0 . 0 
ICOUNT  =0.0 
C 

C  Program  introduces  itself 
C 

WRITE  (*,*)’ 

PRINT  8001 

8001  FORMAT (1X,70( ’*' )) 

PRINT  8002 

8002  FORMAT  (IX , ’ * ’ , 68X , ’ * ’ ) 

PRINT  8003 

8003  FORMAT  (IX,’*  'PCDFRAC'  -  A  PROGRAM  WHICH  CALCULATES  WAVE’ 

1  , IX, ’AMPLIFICATION  FACTORS  *’) 

PRINT  8004 

8004  FORMAT  IX,  ’  »  BASED  ON  THE  COMBINED  DIFFRACTION  AND  REFLECTION’ 

1  , IX, ’ OF  WAVES  BY  A  *’ ) 

PRINT  8005 

8005  FORMAT  (IX,’*  VERTICAL  WEDGE  OF  ARBITRARY  WEDGE  ANGLE.  THE’ 

1  , IX, 'PROGRAM  ALSO  CALCU-  *’) 

PRINT  8006 

8006  FORMAT  (IX,’*  LATES  THE  FHASE  OF  THE  AMPLIFIED  WAVE.  THE  PROGRAM’ 
1  , IX, 'WILL  DEFAULT  TO  * ’ ) 

PRINT  3007 

8007  FORMAT  (IX,’*  ZERO  WEDGE  ANGLE  ( SEMI -INFINITE  BREAKWATER  CASE)’ 

1  .  iX,  No  WEDGE  ANGLE  *’  ) 

PRINT  8008 

8008  FORMAT  (IX. ’+  IS  SPECIFIED .'. 54X ,’*’ 1 
PRINT  8002 

PRINT  8001 
WR I T  E  (+.♦■)  ' 

0 

C  Set  option  for  writing  output  into  file 
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c 

991  PRINT  2199 
2199  FORMAT  <  36 (  '  -*  )  ) 

PRINT  * , ’ DO  YOU  WANT  OUTPUT  IN  A  FILE  ?  (0  =  ’’NO",  1="YES“)  =  ’ 

READ  *,  MM 

IF  (MM  .NE.  1  .AND.  MM  .NE.  0)  GOTO  991 
IF  (MM  . EQ.  0)  GOTO  993 
PRINT  2199 

614  PRINT  *, ’ENTER  NAME  OF  OUTPUT  FILE  -  ’ 

READ  (*, ’ (A) ’ )  INFILE 

INQUIRE  (FILE=INFIL£,EXIST=EXST) 

IF  (EX3T  . EQV .  .TRUE.)  THEN 

613  PRINT  *, ’FILE  ALREADY  EXISTS;  SELECT  ONE  OF  THE  FOLLOWING:’ 

PRINT  *, 1 
PRINT  615 

615  FORMAT (IX, ’0=ENTER  ANOTHER  FILE  NAME  ;  l=OVERWRITE  OLD’ 

1  , IX, ’FILE  WITH  NEW  DATA’) 

PRINT  * , ’ 

PRINT  *, ’OPTION  ?  =  ’ 

READ  *,ANS 

IF  (ANS  .NE.  0  .AND.  ANS  . NE .  1)  GOTO  613 
IF  (ANS  . EQ.  1)  THEN 

OPEN  ( 2 , FILE= INFILE , STATUS^ ’ UNKNOWN ’ ) 

ELSE 
GOTO  614 
END  IF 
END  IF 

OPEN  (  2 , FILE  =  INFILE , STATUS= ’ NEW ’ ) 

WRITE  (2,2201)  INFILE 

2201  FORMAT  (IX,’ OUTPUT  OF  PROGRAM  "PCWEDGE"  -  FILENAME  IS  \Ae) 

993  IF  (ICOUNT  . NE .  0.00)  GOTO  998 

C 

C  Set  Option  for  Interactive  or  Batch  operation 

C 

PRINT  2199 

PRINT  *,’IS  INPUT  DATA  LOCATED  IN  A  FILE  ?  ( 0= "NO" , 1 = " YES" )  =  ’ 
READ  *,  NN 

IF  (NN  .NE.  0  .AND.  NN  . NE .  1)  GOTO  993 
IF  (NN  .EQ.  0)  GOTO  600 
C 

C  Batch  Operation  Option 

C 

PRINT*, ’ENTER  NAME  OF  INPUT  FILE  =  ’ 

READ  ( * , ’ (A) ’ )OUTFILE 

OPEN  i,  i  ,  FILE=OUTFILE ,  STATUS-  uLu'  ) 

998  KCOUNT=0 . 00 
PRINT  2199 
WRITE  (2,2199) 

PRINT  2143. OUTF1LE 
WRITE  (2.2143)  OUTFILE 
2143  FORMAT ( IX, ’ INPUT  FILE  =  ’,A8) 

C 

C  Read  the  file  -  can  have  multiple  values  of  h.t.wavea,  and  wedgea 

C 

K=1 
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1001  READ  ( 1 , 2100,END=999)  D , T , WAVEA , WEDGEA 
2100  FORMAT (4F10 . 2) 

IF  (D  . LT.  0.00)  THEN 

PRINT*, ■WATER  DEPTH  CANNOT  BE  NEGATIVE.  EXECUTION  STOPPED.’ 
PRINT*, ’ PLEASE  CORRECT  INPUT  FILE  AND  RE  -  RUN.' 

GOTO  700 
ELSE 
END  IF 

IF  (T  .LT.  0.00)  THEN 

PRINT* , ' WAVE  PERIOD  CANNOT  BE  NEGATIVE.  EXECUTION  STOPPED.’ 
PRINT*, 'PLEASE  CORRECT  INPUT  FILE  AND  RE  -  RUN.’ 

GOTO  700 
ELSE 
END  IF 

IF  (WEDGEA  .LT.  0.00  .OR.  WEDGEA  .GT.  180.00)  THEN 

PRINT  * , ’ WEDGE  ANGLE  MUST  BE  BETWEEN  0.0  AND  180.0  DEGREES.' 

PRINT  *, ’EXECUTION  STOPPED.’ 

PRINT  *, 'PLEASE  CORRECT  INPUT  FILE  AND  RE  -  RUN.’ 

GOTO  700 
ELSE 
END  IF 

IF  ( KCOUNT  .NE.  0.00)  GO  TO  6551 
WRITE  (2,2199) 

PRINT  2199 
PRINT  2049 
WRITE  (2,2049) 

2049  FORMAT (’  ’) 

WRITE( 2,2050) 

PRINT  2050 

2050  FORMAT (3X, ’DEPTH’ ,3X, ’PERIOD’ ,2X, 'WAV. ANGLE’ ,1X, ’ WAVLNG . ’ , 

1  ’WDG.ANG. ’ ,3X, ’X’ ,7X, ’Y’ ,5X, ’AMP’ ,4X, ’PHA’ ) 

WRI TE ( 2,2051) 

PRINT  2051 

2051  FORMAT ( 4X, ’ (FT) ’ ,4X, ’ (SEC) ’ ,4X, ’ (DEG) ’ , 4X, ’ (FT) ’ ,3X, ’ (DEG) ’ ,4X, 
1  ’ (FT) ’ , 4X, ’ (FT) ’ , 9X, ’ (RAD) ’ /  36 ( ’  -’)/) 

C 

C  Call  subroutine  PADES  to  calculate  wavelength 

C 

6551  CALL  PADES ( L , D , T ) 

WRITE( 2,2150)  D , T , WAVEA , L , WEDGEA 
2150  FORMAT  )  F8 . 2 . 2 X ,  F6 . 2 , 4X ,  F6 . 2 , 2X ,  F7 . 2  ,  IX ,  F6 . 2  ) 

PRINT  2 1 50 , D , T , WAVEA , L , WEDGEA 
J=1 

98  READ  ( 1,2200, END=99)  X ( J ) , Y( J ) . END 
2200  FORMAT (2F10. 2, 25X, 15) 

C 

C  Call  subroutine  MAIN  to  solve  problem  and  deliver 

C  solution 

C 

CALL  MAIN  (F,FABS( J) ,FPHa( J) ,X(J) ,Y( J) , L, WEDGEA, WAVEA) 

WRITE ( 2,2250)  X(J) , Y( J ) , FABS ( J ) ,FPHA(J) 

2250  FORMAT (42X, 2F8 . 2 , 2X , F4 . 2 , 2X , F5 . 2 ) 

PRINT  2250 ,X( J ) , Y( J ) , FABS ( J ) , FPHA( J ) 

X( J i -0 . 00 
Y(.T  )  -0 . 00 
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99 


999 


222 

C 

C 

c 

994 


C 

c 

c 

600 


15 


702 


C 

c 

c 


FABS( J 1=0.00 

FPHA( J ) =0 . 00 

IF  (END  . EQ.  1)  GOTO  99 

J=J  +  1 

GOTO  98 

CONTINUE 

KCOUNT  =  KCOUNT  +1.00 
K=K+ 1 
GOTO  1001 
CONTINUE 

CLOSE ( 1 , STATUS^ ’ KEEP' ) 

CLOSE ( 2 , STATUS= ’ KEEP ' ) 

PRINT  *, ’ 

PRINT  2199 

IF  (MM  .  EQ.  0)  GOTO  994 
PRINT  222, INFILS 

FORMAT ( IX, 'YOUR  OUTPUT  IS  IN  ’,A8) 


Option  to  process  additional  input  files 
PRINT  *, ' 

PRINT  * ,  ’  DO  YOU  HAVE  ANOTHER  INPUT  FILE  ?  ( 0='"NO"  ,  1  =  "YES"  )  = 
READ  *,II 

IF  (II  .NE.  1  .AND.  II  . NE .  0)  GOTO  994 

IF  (II  .EQ.  0)  GOTO  700 

PRINT  * , ' ENTER  NEW  INPUT  FILE  NAME  =  ’ 

READ  ( *,  * (A) ’ )  OUTFILE 

OPEN  ( 1 , FILE=OUTFILE,  STATUS^ ’ OLD ’  ) 

ICOUNT  =  ICOUNT+ 1 . 0 
GOTO  991 


Interactive  option 
PRINT  2199 

PRINT  *, ’OPTION  TO  SET  WEDGE  ANGLE  (DEFAULT  VALUE  =0.00  DEG.)' 
PRINT  2049 

PRINT  * ,  ' SET  WEDGE  ANGLE  ?  ( 0= "NO" , 1=  ’  YES" )  =  ’’ 

READ  *,  LOGIC 

IF  (LOGIC  .NE.  1  .AND.  LOGIC  .NE.  0)  GOTO  15 

IF  (LOGIC  .EQ.  1)  THEN 

PRINT  *, 'ENTER  WEDGE  ANGLE  (DEG.)  =  ' 

READ  * ,  WEDGEA 

IF  (WEDGEA  . GT .  180.0  .OR.  WEDGEA  .LT.  0.00)  THEN 

PRINT  ••.’WEDGE  ANGLE  MUST  EE  BETWEEN  0.0  AND  180.0  DEGREES.’ 

PRINT  2049 

GOTO  702 

ELSE 

END  IF 

ELSE 

WEDGEA  i  0.0 
END  IF 

Get  the  facts:  Initial  values  input  in  this  step 
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PRINT  2199 


PRINT  * , ' SET  INITIAL  VALUES’ 

PRINT  2049 

602  PRINT  *,  'INPUT  WATER  DEPTH  (.FT.)  =  ' 

READ  * , D 

IF  (D  . LT.  0.00)  THEN 

PRINT  *.  ’WATER  DEPTH  CANNOT  BE  NEGATIVE.  ' 

PRINT  2049 
GOTO  602 
ELSE 
END  IF 

603  PRINT  'INPUT  INCIDENT  WAVE  PERIOD  (SEC.)  -  ’ 

READ  *,T 

IF  (T  .LT.  0.00)  THEN 

PRINT  *, ’INCIDENT  WAVE  PERIOD  CANNOT  BE  NEGATIVE’ 

PRINT  2049 
GOTO  603 
ELSE 
END  IF 

701  PRINT  *,’ INPUT  INCIDENT  WAVE  ANGLE' 

PRINT  * , ’ ( DEG .  CCW  FROM  POSITIVE  X  -  AXIS)  =  ’ 

READ  * , WAVEA 
PRINT  2199 
WRITE  (2,2199) 

C 

C  Call  subroutine  PADES  to  calculate  wavelength 
C 

730  CALL  PADES ( L, D, T) 

IF  (WEDGEA  . NE .  0.00)  THEN 
WRITE  (2,2001)  WEDGEA 

2001  FORMAT  ( IX ,’ OBSTACLE  IS  A  VERTICAL  WEDGE  OF  ANGLE  ’ , F6 . 2 
1  ,  ’  DEGREES '  ) 

ELSE 

WRITE  (2,2002) 

2002  FORMATt IX,  ’ SEMI  -  INFINITE  BREAKWATER  (WEDGE  ANGLE' 

1  , IX, ’ =  0.00  DEGREES) ’ ) 

END  IF 
PRINT  2049 
WRITE  (2,2049) 

WRITE  (2,2050) 

WRITE  (2,2051) 

WRITE  (2,2150)  D , T , WAVEA , L , WEDGEA 
710  PRINT  * , ’ HOW  MANY  POINTS  DO  YOU  WISH  TO  ENTER  (UP  TO  99)?  =  '  • 

READ  *■ ,  KK 

PRINT  667 

667  FORMAT ( IX, 'AT  EACH  PROMPT,  DO  THE  FOLLOWING  :’,/ 

1  , 4X , ' 1 )  TYPE  IN  X  -  COORDINATE’ ,/,4X, '2)  TYPE  A  COMMA' 

2  ,  /  ,  4X ,  ’  3  )  TYPE  IN  Y  -  COORDINATE ’,/,  4X  ,’ 4  )  PRESS  "ENTER’,'/) 

DO  10  1=1, KK 
PRINT  2UOU , 1 

2000  FORMAT  (1X,’(X,Y)  COORDINATES  FOR  POINT  NO.  ’,12 
1  . IX, ’  (FT. )  = 

READ  * . X  (  I )  ,  Y  (  I ) 

C 

C  Call  subroutine  MAIN  to  solve  problem  and  deliver 

0  solution 


1 
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n  o  n 


i 


C 

CALL  MAIN  t K , FABS( 1  )  , FPHAl i > ,  v i i ) . Y ( I ) , L , WEDGEA , WAVEA ! 

WRITE  !  2.2250)  X(1  )  ,  Y(  I  )  .  FABiS >  I  »  ,  FFHA(  I  j 
10  CONTINUE 

CLOSE ( 2 , STATUS- ’ KEEP 1  ) 

FRINT  2199 

IF  ( WEDGEA  NE.  U . 00  )  THEN 
PRINT  LOO  1 , WEDGEA 
FRINT  2049 
ELSE 

PRINT  2U02 
PRINT  2049 
END  IF 
PRINT  2050 
PRINT  2051 

PRINT  2 1 50, D,T, WAVEA. L, WEDGEA 
Screen  Output 
DO  20  N -  1 , KK 

PRINT  2250, XlN) , YIN) , FABS( N ) , FPHAl N ) 

20  CONTINUE 

PRINT  * . ' 

IF  (MM  . EQ.  0)  GOTO  700 
PRINT  222, INFILE 
700  STOP 

END 

C 

(1  Suhr™n  i  [)<■  FADES,  to  calculate  wavelength  . 

C  Based  _>n  Fade’s  approximation  soiution  for  the  linear 

C  dispersion  relation 

C 

SUBROUTINE  FADES ( L , D , T ) 

REAL  L 

PI  -3. 14159 

PL2-2 . *  PI 

..■MG  A -i' 12/ 7 

OMGG  ' :  jMG  A  * *■  2  .  /  ,1 2 . 2 

A  ‘-OMGG-*  D 

BB-A* ( A>1 . / i 1 . t A*( 0 . 66667  +  A* ( 0 . 35550 +  A+( 0 , 16084  + AM  0 . 06320 ♦ A* 

1  (1 . 02174  *AM  0 . 00654  +  A*  (0 . 00171  ♦  A*(0 . 00039+A*0 . 00011  ))))))))  I  i 

b  -  sort i pb i 

PK-B/Ij 
;  -:PI i'K 
RETURN 
END 
c 

i  ‘  Ltibr.'uti  ne  MAIN,  which  normalises  the  input  and  '-ills 

U  ,,utr  ut  i  ne  WEDGES,  which  solves  the  actual  prop  lent . 

ij 

SUBROUTINE  MAINi  F.  FABS  .  FFHA  .  X  ,  V  ,  L ,  WEDGEA  .  WAVEA  ) 

REAL 

PI .7. 14159 

I  I  2  .  r  [-  I 

XX  A ■  l  i 


YY - ABS i Y/L ) 

XRL = SORT ( ( XX*  *2 . 0  )  +  ( YY**2 . 0 ) ) 

IF  iX  .  c.w.  .• .  uO  AND.  Y  .  LT .  O.UU)  THEN 
XTH=1. 5‘FI 
GOTO  351 
ELSE 

IF  (X  .  EQ.  0.00  .AND.  Y  . GT .  u.OO)  THEN 
XTH  =  0 . 5  *  P I 
GOTO  351 
ELSE 

GOTO  4 LS 
END  IF 
END  IF 

428  IF  (XX  .EQ.  u.OO  .AND.  YY  .EQ.  0.00)  THEN 
XTH  -  0.00 
GOTO  351 
ELSE 

XTH=ATAN ( YY/XX ) 

END  IF 

IF  a  .LT.  U.OOIGOTO  100 
IF  ( Y  . LT.  0.00)  GOTO  300 
GOTO  400 

100  IF  (Y  .LT.  0.00)  GOTO  200 
XTH=PI -XTH 
GOTO  400 

200  XTH  =  P I +XTH 

GOTO  400 

30 0  XTH=PI 2-XTH 

400  IDX-0.0 

351  CALL  FCWEDGES ( F , FABS , FPHA , XRL , XTH , WEDGEA , WAVEA, IDX ) 

RETURN 

END 

SUBROUT  I NE  PCWEDGES ( F . FABS , FPHA , XRL , XTH , WEDGEA , WAVEA . 1 DX ) 

0  *  ****  **********  ******■***,**  *  +  +  *****♦ 
C  «  THIS  COMPUTER  PROGRAM  WAS  PREPARED  UNDER  THE  EFFORT  OF  CIVIL  WORK 
C  *  R&D  PROGRAM  OF  COE.  NEITHER  ANY  OF  AGENCIES  NOR  amvqnE  ASSUMES 
’  *  AN i  -SGAL  LIABILITY  OR  RESPONSIBILITY  FOR  THE  ACCURACY  OF  THE 

*  PROGRAM. 

C  *  +  **•****.*..  *<:*********  +  *»****.*•*♦♦♦ 


WAVE  SCATTERING  BY  A  WEDGE  OBSTACLE. 
C  NN.  IDX-u  TO  CALL  BESJ . 


NUMBER  OF  SUMMATION  TERMS  IS 


I  •iP.AM7.TER  NN  -  LOO 
b I  MENS  IN  BJ  (  NN  i  ,  W(NN  )  .XM(NN) 
COMPLEX  F.TM 

LATA  Is LR/ 1 . E-8/ , ITER/8/ 

IF  <  I  LX  NE.  •>.  GOTO  4 
Lu  ;  .  -  I  ,  NN 

L'  J  !  I  .1  -  > .  .  I  lO 

■•.nt:n"e 

:  i  -  s  .  I  4 1 5  jL 8 54 
CPI -P i / 180  . 

XKR-XRL»L  ..‘PI 
IH-XI!; 


B7 


WA-WAVEA+CFI 

XNU  =  (  J60  .  -  WEDGEA  i  180. 

XH (ii  1  .  u 

IF  ( 1  LX . HE . 0 i  GOTO  14 
CALL  BE5J  4  XKR ,  U  .  0  ,  1 ,  W  ,  NZ  ) 

B  J  (  1  i  =  W  (  1  ) 

I COUNT -0 
DO  11)  N  =  1  ,  NN 
N1=N+1 

I  COUNT : I COUNT+1 

IF ( I COUNT. LE. ITER)  GOTO  6 

MNN-N 

GO  TO  *  4 

a  XM( N 1 ) =  FLOAT ( N ) /XNU 
M= INT ( XM ( N 1 )  ) 

ALPHA-XMlNl ) -M 
Ml-M+1 

CALL  BESJ( XKR, ALPHA, Ml. W.NZ) 

IF(Nl.EQ.NN)  WRITE ( 6,9)  XKR, ALPHA , Ml , W( Ml ) , NZ 
N  FORMAT)/'  ****  NO.  OF  TERMS  FOR  SUMMATION  IS  INSUFFICIENT  ***■*’, 
1  /’  XKR, ALPHA, Ml ,W(M1 ) ,NZ  a '  , 2F10 . 4 , 15 , E15 . 6 . 11/ ) 

BJ  (  L'l  )  -W(  Ml  ) 

IF;  ABS I  BCi.Nl)).  GT.  TOLRj  I  COUNTS 
10  CONTINUE 
14  CONTINUE 
F  a  B  J  (  1 

DO  f:u  N  a  1 ,  NNN 
N1  CM 
XMN-XM( N1  i 

TM=  i-i.j.  1 . 0  )  **XMN-*BJ(NI  )  ♦COS(XMN*WA ) +COS( XMN*TH ) 

F  -  p  t  TM 

LIU  CONTINUE 
'•  4  XNU  *  F 
r  k -  REAL , F  ) 

FI -AIMAG( F ) 

FABS-COKTi FR*FR+FI*FI ) 

: F : WA . LE . 1 . z- 8 )  FABC -FABS/2 
'  Ft  FAbS  .  L  i' .  TOLR  i  GOTO  30 
1  PHA  -  AT  AN2 i F I .  FR  ) 

BET 'TCI 
3u  BPHA-U  I; 

RETURN 

END 

:’BR.  ''TTNE  !  ECU.  v  .  ALF  HA.  N  .  Y  .  N2  ) 

.'■ANLI A  MATHEMATICAL  PROGRAM  LIBRARY 
APPLIED  MATHEMATICS  DIVISION  2642 
■  'AN I  IA  LABORATORIES 
:-~jd  •«|JiR<y  .  NEW  MEXICO 

'ONTRi >L  DATA  8 A 00  VERSION  8.1  JANUARY  1076 
*  *  *  '  **»*»****♦***♦•*,  ***♦*«♦*  ♦*«**(. 

'  ISSUED  BY  SANDIA  LABORATORIES, 

•  A  PRIME  CONTRACTOR  TO  THE 

*  'Nil  EL  •  -TATES  ENERGY  RESEARCH  AND  DEVELOPMENT  ADMINISTRATION 

NOTICE  *  *  *  *  *  *■  •  *  *  *  *  *  *  * 


APPENDIX  C:  NOTATION 


I 

I 

I 


I 

I 


« 


« 


A 

o 

a 

o 

e 

8 

h 

i 

J 

o 


k 


(r , 0 , z) 
r 
u 
z 
a 
8 
n 

ni 

e 

6 

o 

V 

Tf 


-igaju 

Incident  wave  amplitude 

2.71828....;  base  for  natural  logarithm 

Gravitational  acceleration 

Water  depth 

v^T 

Zeroth  order  Bessel  function  of  the  first  kind 
“  order  Bessel  function  of  the  first  kind 
Wave  number 

Cylindrical  coordinate  system 
Radial  coordinate 
Flow  velocity 
Vertical  coordinate 

Approach  angle  of  incident  wave  train 
Phase  of  <Kr,0) 

Water  surface  elevation 
Incident  water  surface  elevation 
Angular  coordinate 

Angular  coordinate  for  water  domain 

0  /  2tt 
o 

3.14159 . 


4>(r,0,z,t) 

f(r,0) 


Three-dimensional,  time-dependent  velocity  potential 
Velocity  potential  in  horizontal  plane 
Radian  frequency 


Subscripts : 


r 

p 

Mathemat ical 

3 

r 

t  I 

|  j 

Tmf  ) 
Ref  ) 


Denotes  r-coordinate  variable 
Denotes  6 -coordinate  variable 
symbols : 

Partial  differentiation 
Summat ion 
Absolute  value 
Imaginary  part 
Real  part 


Cl 


